The spin anisotropy of the magnetic excitations in the normal and superconducting 
states of optimally doped YBa 2 Cu 3 06.9 studied by polarized neutron spectroscopy 
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We use inelastic neutron scattering with spin polarization analysis to study the magnetic exci- 
tations in the normal and superconducting states of YBa2Cu306.9- Polarization analysis allows us 
to determine the spin polarization of the magnetic excitations and to separate them from phonon 
scattering. In the normal state, we find unambiguous evidence of magnetic excitations over the 
10-60 meV range of the experiment with little polarization dependence to the excitations. In the 
superconducting state, the magnetic response is enhanced near the "resonance energy" and above. 
At lower energies, 10<i5<30 meV, the local susceptibility becomes anisotropic, with the excitations 
polarized along the c-axis being suppressed. We find evidence for a new diffuse anisotropic response 
polarized perpendicular to the c-axis which may carry significant spectral weight. 



I. INTRODUCTION 

High temperature superconductivity (HTS) arises 
when certain two dimensional antifcrromagnctic Mott in- 
sulators are electron or hole doped [1]. The antifcrro- 
magnetic parent compounds such as La2Cu04 show spin- 
wave excitations up to 2J«300 meV [2]. Doping causes 
the magnetic response to evolve from that of spin waves 
to a more structured response [3-13], with strong spin 
fluctuations being observed for superconducting compo- 
sitions in a number of systems including YBa 2 Cu306+:i; 
(YBCO) [7-10, 14, 15], La 2 - K Sr 2 ,Cu0 4 [11, 16] and 
Bi 2 Sr 2 CaCu20 8+( 5 [12, 13, 17]. Many optimally-doped 
cuprates show a strong well-defined collective magnetic 
excitation which is localised in reciprocal space and 
strongest near the Q=(l/2,l/2)=(7r, it) position. It is 
sharp in energy and develops on cooling through the crit- 
ical temperature. This excitation has become known as 
the "magnetic resonance". The magnetic resonance has 
been observed in YBa2Cu306+a: [4-6], Bi 2 Sr 2 CaCu 2 08+5 
[17], Tl 2 Ba 2 Cu0 6+<5 [18] and HgBa 2 Cu0 4 +5 [19]. 

The magnetic resonance is certainly the strongest fea- 
ture in the magnetic excitations spectrum of the mate- 
rials listed above, however, it only accounts for a small 
faction («2%) [9, 15, 17] of the total scattering expected 
from the unpaired 3d electrons of the Cu atoms. In this 
work we search for other contributions to the response 
which are spread out in energy and wavevector but nev- 
ertheless may carry significant spectral weight. These 
are harder to observe because they are weak and may 
not show the strong temperature dependence which al- 
lows the resonance to be easily isolated. We use inelastic 
neutron scattering with polarization analysis to isolate 
the magnetic scattering from phonon scattering. 

We find that there is a significant response in the nor- 
mal state which can account for much of the spectral 
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weight from which the resonance is formed. In the su- 
perconducting state, we find evidence for a diffuse con- 
tribution at energies well below the resonance. This new 
contribution is polarized with strong fluctuations perpen- 
dicular to the c-axis. 



II. BACKGROUND 
A. Polarization Analysis 

Neutrons scatter from condensed matter via two pro- 
cesses: (i) The electromagnetic interaction probes fluc- 
tuations in the magnetization density of the electrons 
(in this paper this is referred to as magnetic scattering), 
(ii) The strong nuclear force is responsible for scattering 
from the atomic nuclei. The nuclear scattering allows 
us to probe phonons which are correlations (in time and 
space) between the position of the nuclei. The existence 
of two distinct scattering processes makes the neutron an 
extremely versatile probe. However, it also means that 
the two types of scattering can mask each other. 

Polarization analysis of the neutron's spin allows the 
separation of magnetic and nuclear (phonon) scattering. 
In the present work, we use longitudinal polarization 
analysis (LPA). In LPA, a spin-polarized incident neu- 
tron beam is created and its polarization maintained by 
a small magnetic field (~1 mT). The number of neu- 
trons scattered with spins parallel or antiparallel to this 
quantizing field are then measured. We label each spin- 
polarization state as parallel (=up,f,+) or antiparallel 
(=down,J,,— ) to the applied field. The cross sections are 
referred to as spin-flip (SF) (f— kU4- — ^t) or non-spin-flip 
(NSF) ('['— >t,4_— A natural reference frame in which 
to understand the cross sections is one referenced to the 
scattering vector Q = k; — ky of the neutron, where 
kj and kj are the incident and final wavevectors of the 
neutron. Thus, x || Q, y _L Q, and i _L Q and _L to 
the spectrometer scattering plane (the plane containing 
kj and k/). We make measurements with the neutrons 



polarized along each of these axes. 

The neutron cross sections as a function of spin polar- 
ization have been derived and presented elsewhere [20- 
24]. The spin- flip magnetic cross section for spin polar- 
ization || Q is 
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where (7r o ) 2 =0.2905 barn sr _1 and |F(Q)| 2 is the 
anisotropic magnetic form factor for a Cu 2+ d x % i or- 
bital. x"^(q ; ^ w ) is the generalized susceptibility cor- 
responding to magnetic fluctuations along the u-axis. 
Thus, for example: 
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where the angle brackets denote thermal averages. The 
spin-dependent cross sections including the nuclear co- 
herent cross sections (i.e. the phonon cross section) 
N(q,w) are: 
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where we have neglected the nuclear spin incoherent 
cross-section which is small in the present experiments 
[25] and BG denotes the background for the configura- 
tion. In this work we isolate two components of the sus- 
ceptibility by comparing different SF cross sections: 
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B. Bilayer Effects 



YBa2Cu306+iE has two Cu02 planes per unit cell (See 
Fig. 1). The usual starting point for models of the mag- 
netic response is to neglect the electronic coupling be- 
tween Cu02 planes in different unit cells and include only 
coupling between the Cu02 planes of the bilayer located 
in the center of the unit cell in Fig. 1. This leads to a 
pair of bonding (b) and antibonding (a) energy bands. 
The presence of a mirror plane between the two planes 
of the bilayer means that the magnetic excitations have 
distinct odd (o) or even (e) character. In this description, 
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FIG. 1. (Color online) The ortho-I structure of YBa2Cus06 



the magnetic response is of the form [26-29] 
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where d is the separation of the Cu02 planes. For 
YBa 2 Cu 3 06.g d=3.38 A, this means the odd re- 
sponse is strongest at l=(n + l/2)c/(2d)=1.73, 5.3,... 
The strongest features in the magnetic response of 
YBa2Cu306+x observed by INS are in the odd channel 
[4-6] and we measure the odd channel in the present ex- 
periment. We note that weaker resonance features have 
been reported in the even channel [30, 31] for various 
dopings. The reported even resonance occurs at higher 
energy than in the odd channel. 



C. Sample Details 

We investigated a near optimally doped sample of 
YBa 2 Cu 3 06.g (T c =93 K) grown by a top seed melt 
growth technique [32]. YBa 2 Cu 3 06.9 has the ortho- 
I structure show in Fig. 1 with lattice parameters 
a=3.82 A, 6=3.89 A and c=11.68 A (T=77 K) [33]. 
The single crystal studied in the present experiment is 
twinned and the results presented are an average over 
the two twin domains. The crystal had a mass of 32.5 g 
and mosaic spread 1.3°. It was annealed for 17 days 
at 550°C, followed by 13 days at 525°C, in oxygen to 
achieve the required oxygen stoichiometry. Neutron de- 
polarization measurements (see Fig. 4) indicated that 
T c (onset) = 93±0.2 K. Based on T c and the heat treat- 
ment [34, 35], we estimate the oxygen stoichiometry to 
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FIG. 2. (Color online) (a) Illustration of the reference frame 
used to describe polarization analysis. Neutrons are po- 
larized along the x, y or z axes. x is parallel to Q, 
z is perpendicular to x and in the (h, k, 0) plane. Thus 
x || (1.5a*+0.5b*+1.8c*), y || (0.54a*+0.18b*-4.6c*), z || 
(-0.5a* + 1.5b*), 9 = 20.6° and cos 2 9 = 0.88. (b) Illustra- 
tion of the area in reciprocal space where the measurements 
in Sec. IIIB were made. For E < 52 meV, we used data 
collected over the black square (1.3</i<1.5 and 0.3<fc<0.5) 
to infer x"^) measured over the grey area. Data in Fig. 6 
covers the black square plus dotted area. 



be x=0.9±0.01. 



D. Experimental Setup 

Experiments were performed using the IN20 three-axis 
spectrometer at the Institut Laue-Langevin, Grenoble us- 
ing a standard longitudinal polarization analysis set up. 
Neutron polarization analysis was carried out using a fo- 
cussing Heuslcr monochromator and analyzer. The sam- 
ple was mounted with the [310] and [001] directions in the 
horizontal scattering plane of the instrument. We worked 
around the (1.5,0.5,1.73) reciprocal space position so as 
to avoid strong phonon scattering near E «40 meV [6]. 
We used a sample goniometer to access reciprocal space 
positions out of the (3/i, h, I) plane. Data were con- 
verted to an absolute scale using a vanadium standard 
and Eq. 1 and measurement of an acoustic phonon at 
Q=(0.2,0.2,6). The overall error in the absolute scale is 
about 20%. We use the reciprocal space of the average 
tetragonal lattice (with a ps 3.86 A) to label wavevectors 
with Q = ha* + kb* + Ic*. 

In order to reduce neutron depolarization for measure- 
ments made in the superconducting state, the sample 
was cooled through T c and to 10 K while shielded by 
a /U-metal shield such that (j-qH <0.3 fiT. During the 



measurement, fields in the range \^qH.\ =0.7-0.11 mT <c 
/io-ffci~25-85 mT were applied to the sample. Therefore, 
the sample was in the Meissner state. 

The finite polarization of the incident neutron beam 
and other instrumental imperfections leads to a mixing 
of the spin-flip and non spin-flip channels. This can be 
described by a flipping ratio F, where the measured cross 
section is: 
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We corrected our data for this mixing using the standard 
equations [36]: 



coir 
°SF 



F 



F - 1 



'SF 



1 



F - 1 



Trie as fi-r\ 
CT NSF ) {<) 



where the flipping ratio F « 7.5 was determined from 
measurements on Bragg peaks made under the same con- 
ditions. For experimental reasons, measurements were 
made with neutrons polarized parallel and perpendicular 
to the scattering vector Q which meant that the neutron 
polarizations and hence the measured susceptibilities are 
not along the crystallographic axes (see Fig. 2). For ex- 
ample, the angle between the y-axis and the crystallo- 
graphic c-axis is 0=20.6°. This leads to a small mixing 
of the different components of the susceptibility during 
the measurement. Thus: 

4i - K QA1 Xa + 0-01X6 + 0.88 X ^'(= x") 

ali - all * O.lx" + 0.9x1; (= x"/b)- (8) 

This mixing does not affect the conclusions of the paper 
and we have not corrected for it. We refer to the two 
components above as x'a/b and x" c . The local susceptibil- 
ity (see Sec A 1) was estimated by measuring a grid of 36 
points over the area 1.3 < h < 1.5 and 0.3 < k < 0.5 
(at the highest energy we used 1.5 < h < 1.7 and 
0.5 < k < 0.7 in order to close the scattering triangle). 
Points were weighted according to the number of equiv- 
alent positions in the grey area of Fig. 2(b). 



III. RESULTS 

A. Energy- and Wavevector-Dependent Scans 

Fig. 3 shows energy-dependent scans made at the 
(1.5,0.5,1.73) position with various spin polarizations. 
At this position in reciprocal space the non spin-flip 
(phonon) scattering is up to 8 times stronger than the 
spin-flip scattering. Thus an unpolarized measurement 
made under the same conditions would be dominated 
by phonon scattering at some energies (the comparison 
with unpolarized experiments is discussed further in Ap- 
pendix B). In the normal state the <r|* cross section 
is larger than a^ and crj* over a wide energy range, 
20 < E < 60 meV, signalling the presence of magnetic 
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FIG. 3. (Color online) Energy-dependent scans with polar- 
ization analysis collected at Q=(1.5,0.5,1.73). (a-b) Spin-flip 
and non-spin-slip cross sections for various spin polarizations 
in the superconducting (T = 10 K) and normal (T = 94 K 
— T c + 1 K) states, (c) Out-of-plane generalized suscepti- 
bility Xc determined from (a) and (b). The solid line for the 
T = 10 K data is a resolution corrected fit to the cross-section 
described in the text. The horizontal bar represents the full- 
width-at-half maximum (FWHM) resolution for a S(oj — u>o) 
cross section. 



excitations. We can use Eq. 8 to isolate the out-of-plane 
response x'li this is shown in Fig. 3(c). In the super- 
conducting state there is a large increase in erjj: and (7^4 
(cr|| was not measured in this case) near the resonance 
energy. The difference scan Fig. 3(c) shows a sharp res- 
onance peak at E 40 meV which appears to have 
formed by a transfer of spectral weight from lower ener- 
gies E < 35 meV. The x" response appears to be largely 
gapped below about 30 meV. Similar data was obtained 
using unpolarized neutrons by Bourges et al. [37]. We 
do not observe a collective magnetic excitation in the 50- 
60 meV range as observed recently in HgBa 2 Cu04 +l 5 [38]. 
We note that there is a peak in the non spin flip channel 
in this energy range in Fig. 3(a). 

In order to analyze our data further, we fitted the T = 
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FIG. 4. (Color online) (a)-(f) Wavevector-dependent scans 
with LPA through the (1.5,0.5,1-73) position at the resonance 
energy (a,b) and lower energies (c-f). The solid lines are reso- 
lution corrected fits to Eq. 10. The horizontal bar represents 
the FWHM resolution for a 5(Q — Qo) cross-section. The in- 
set to (b) shows T c determined by a neutron depolarization 
technique in which the sample was field cooled through T c in 
a vertical field. The field was then rotated to be horizontal 
and the spin flip scattering on the (310) Bragg peak measured 
on warming. 



10 K scan in Fig. 3(c) to the resolution-corrected model 
cross section 



Xc(q, w) = [A5(uj - u ) + B6(uj - oj q )} 
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where 6 is the heaviside step function and a is the width 
parameter extracted from a q-dependent scan through 
the resonance (see Table I). Throughout this paper we 
use the RESTRAX simulation package [39] to perform 
convolutions of the instrumental resolution function and 
model cross sections. Using the cross section defined by 
Eq. 9, we find that the width of the peak due to the 
resonance in Fig. 3(c) is resolution limited and Hloq 
41 ■ 1 meV. 

We have converted the data in Fig. 3(c) to absolute 
units using Eq. 1 without attempting to deconvolve the 
experimental resolution. This means that each point in 
the scan is an average (in wavevector and energy) over 
the instrumental resolution. Keeping this in mind, we 
have integrated the response in Fig. 3(c) in energy for 
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FIG. 5. The out-of-plane response x" in the normal and su- 
perconducting states determined from the data in Fig. 4. The 
solid lines are resolution corrected fits to the cross-section de- 
scribed in the text. 

4 < E < 60 mcV for T=10 K and 94 K. From Eqs. 2 and 
A5, we find the out-of-plane fluctuating moments (m^.) 
are 0.50±0.05 and 0.48±0.05 ^lur 1 at T=10 K and 
94 K respectively (these are averaged over the resolution 
width in wavevector shown in Fig. 4). Thus this increase 
in the response at the resonance energy can be accounted 
for by a shift in spectral weight from lower energies. 

Fig. 4 shows wavevector dependent scans along the 
(1.5, k, 1.73) line at three characteristic energies. Fig. 5 
shows the susceptibility extracted from the data in Fig. 4 
using Eq. 4. In the normal state (T = 94 K), we observe 
a magnetic response at all three energies. On cooling to 
T = 10 K, the lower frequency E = 26 meV response 
is suppressed while the response at the resonance energy 
(E = 40 meV) increases dramatically and the q-width de- 
creases. The data were fitted to a model consisting of four 
incommensurate peaks with locations Qs = (l/2±5, 1/2) 
and (1/2, 1/2 ± 6) and width a: 

x"(q^)=A^exp{- (Q 2 y 2 }. (10) 

The results of this fitting procedure are shown in Table I 

We first consider the scans at the resonance energy 
(hu=40 mcV). A single Gaussian peak (<5=0) provides a 
good description of the scan in the superconducting state 
[Fig. 4(a) and Fig. 5(a)]. In the normal state, there is 
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TABLE I. Incommensurability 5 and width a parameters 
obtained from fitting Eq. 10 to the scans in Fig. 4. Where no 
error is quoted, the parameter was fixed. 



magnetic scattering at the resonance energy [Fig. 5(b)]. 
The existence of a magnetic response at this energy in 
optimally doped YBCO has been a subject of some de- 
bate [4-6, 37, 40] and we will discuss this later. It is clear 
from our data that the response at the resonance energy 
is broader in q and weaker in the normal state than the 
superconducting state. If we fit the 40 mcV data using 
Eq. 10 with 6=0, wc find cr=0.18±0.02 and 0.115±0.01 
for the normal and superconducting states respectively. 
Returning to the superconducting state data at lower en- 
ergy, we find a single Gaussian peak (<5=0) does not pro- 
vide a good description of the ftw=34 meV (T=10 K) 
scans [Fig. 4(c) and Fig. 5(c)] in the superconducting 
state. Better fits are obtained when a finite incommen- 
surability <5=0.12±0.02 is used. This i5 is in agreement 
with that obtained in other studies of optimally doped 
YBCO [15, 40]. In the normal state [Fig. 4(b,d,f) and 
Fig. 5(b,d,f)] we see clear magnetic scattering at the three 
energies investigated. We do not see a two-peaked struc- 
ture as in Fig. 4(c), instead the response appears to be 
broadened out into single peak which, in some cases [e.g. 
Fig. 4(b,f)], looks "fiat topped". To contrast the nor- 
mal and superconducting state responses, we have fit- 
ted the scans with the value of 6 determined from the 
T=10 K and ftw=34 meV scan. The normal state re- 
sponse is broader in all cases (see Table I). 



B. Local Susceptibility Measurements 

In order to search for the diffuse contributions to the 
magnetic response, we sampled a grid of points near 
the (3/2,1/2) position where the response is generally 
stronger. Extended grids at two characteristic energies 
are shown in Fig. 6. For this part of the experiment 
we collected three spin-flip channels and we were able 
to extract x'a/t anc ^ Xc- The l° west row of Fig. 6 shows 
the signal extracted via Eq. 4. The data collected at 
E = 40 meV shows that the response is strongest near 
the (1.5,0.5,1.73) position both in the normal and su- 
perconducting states. At E = 26 meV, we see a nor- 
mal state response which is spread out: see, for example, 
x" /b (E = 26 meV,T = 100 K), where the upper part of 
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FIG. 6. (Color online) (Top two rows) Spin polarized cross sections as defined in the text for E=26, 40 meV and T—10, 100 K. 
The tt channel is mostly phonon scattering and the f4- channels are mostly magnetic scattering. (Bottom row) Magnetic 
responses Xa/bil' w ) ano - Xc( t l: £ *') determined from top two rows. Note that some structure is due to statistical noise. 
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FIG. 7. (Color online) (a)-(b) Integral of the spin-polarized 
cross sections over reciprocal space described in text and 
Fig. 2. (c)-(d) Local susceptibility determined from data in 
(a)-(b) via Eqs. 4 and Al. 

the map shows signal. On entering the superconducting 
state x'c shows a much larger change than x'a/b suggesting 
that a spin anisotropy develops. 



Fig. 7 shows the wavevector integrals collected at a 
number of energies over the grey region shown in Fig. 2. 
This is the region of highest intensity in the Brillouin 
zone, but there is clearly scattering in other parts of the 
zone. The contribution of the grey region to x'oi 1 ^) is 
shown in Fig. 7(c) and (d). Fig. 7 shows that there is a 
strong response in the normal state over a wide energy 
range. When compared to the energy- dependent scan 
at (1.5,0.5,1.73), we see that the higher energy response 
is relatively stronger. This is due to the presence of a 
broader response in q at higher energies E > 50 meV 
[10, 15, 41]. On entering the superconducting state, we 
see a strong reduction in x'c with little change in x'a/b- 
This means the magnetic response develops a strong spin 
anisotropy in the superconducting state (see Sec. IV B 
for more discussion). For higher energies, E > 34 meV, 
the response increases in the superconducting state, not 
only at the resonance energy, but up to 60 meV. Table II 
shows that when integrated over the range 12 < E < 
60 meV the total fluctuating moment (m 2 ) increases by 
about 60%. In order to compare with other studies of the 
resonance in near optimally doped YBCO [9, 15, 17], we 
have also integrated the data in Fig. 7 over the smaller 
energy range 30 < E < 60 meV (see Table II) in this case 
we see a larger change in (m 2 ) (between the normal and 
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12 < ftw(meV) < 60 




10 


0.031 ± 0.004 0.026 ± 0.004 


0.088 ± 0.007 


100 


0.017 ± 0.003 0.022 ± 0.003 


0.056 ± 0.005 
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0.024 ±0.003 0.026 ±0.003 
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100 
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TABLE II. Fluctuating moments {m 2 a / b ), (m 2 ,) and (m 2 ) = 

2( m a/b) + ( m c) m the normal (T = 100 K) and superconduct- 
ing (T = 10 K) states calculated by numerically integrating 
the response in Fig. 7. The errors quoted are statistical and 
do not include the systematic error in the absolute scale which 
is about ±20%. 
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FIG. 8. Schematic representation of \' (c\,uj) in the supercon- 
ducting state of YBa2Cu306.9 based on Refs. 42 and 43. The 
black line is the resonance mode and grey area the particle- 
hole continuum. Scans (a), (b) and (c) correspond approxi- 
mately to 20, 40 and 60 meV. 



superconducting states) which is comparable to previous 
reports [9, 15, 17]. 



IV. DISCUSSION 

A. Response in the Normal and Superconducting 
States 

Theories of the magnetic excitations in the super- 
conducting state of cuprate superconductors such as 
YBa 2 Cu 3 06+x are well developed [42-53]. Many fea- 
tures are explained by a magnetic exciton scenario 
[42, 43, 46, 47, 52] in which the resonance is a bound 
state in the particle-hole channel, which appears in a re- 
gion of q — u/ space where there are no damping processes 
due to electron-hole pair creation. This is illustrated 
schematically in Fig. 8. In such a picture, significant 
magnetic response should also be present in the normal 
state. As the system enters the superconducting state 
we expect the low energy response to be suppressed be- 
low E < A and an enhancement of the response at the 



resonance energy. This is the behaviour seen in Figs. 3 
and 5. The nature of the magnetic response in the nor- 
mal state of optimally doped YBCO has been a subject 
of debate, particularly with regard to energies near the 
resonance energy [4-6, 37, 40, 54]. Some studies suggest 
there is a significant response [4, 37] for q w (1/2, 1/2) 
and fko ~ 40 meV, while others claim the response is 
absent or too weak to observe [6, 40, 54]. The present 
experiment allows the magnetic response to be separated 
from phonon scattering. We find that the out-of-plane 
response ^'(qjw) is peaked around frw w 30 meV for 
q « (1/2, 1/2) in the normal state (T = 94 K). On cool- 
ing there is a shift of spectral weight to higher energies 
which leads to the formation of the resonance peak near 
40 meV, with the concomitant formation of incommen- 
surate peaks observed at 34 meV and a spin gap below 
about 30 meV for the x" component of the response. This 
is consistent with the formation of a magnetic excitonic 
mode as illustrated schematically in Fig. 8. The work 
presented in this paper refers to optimally doped YBCO 
where it is harder to separate the magnetic contribution 
from phonons and other background scattering than for 
other compositions. We note that for undcrdoped YBCO 
(e.g. YBa 2 Cu306.e) [8, 10, 55, 56] a strong dispersive 
excitonic mode is also observed in the superconducting 
state. On warming to T c the remnants of this mode are 
clearly observable and persist well above T c . 

The discussion above relates to the energy- and 
wavevector- dependent scans presented in Sec. Ill A. 
These yield information about the out-of-plane fluctua- 
tions described by x" ■ We did not collect the correspond- 
ing scans for x'a/b' however, we did probe this compo- 
nent of the local susceptibility in the measurements pre- 
sented in Sec. IIIB. These measurements were designed 
to yield estimates for the total response in a region of q 
space rather than identify the location of specific features 
such as incommensurate peaks. They are summarized 
in Fig. 7(c) and (d). In Fig. 7(c) we see that there is 
strong evidence for additional scattering below 30 meV 
in the x'a/b component of the response. This response 
appears to be rather spread out in wavevector when we 
inspect the corresponding map (huj =26 meV, T — 10 K) 
in Fig. 6. Thus our results suggest that there are other 
(diffuse) contributions to the x'a/b response at low ener- 
gies in the superconducting state. The x'a/b component 
of the response has a lower 'spin gap' than the x'c com- 
ponent. The low energy response (E < 30 meV) may be 
due to the electron-hole continuum also present in the- 
ories of the resonance [42, 43, 52]. This is illustrated 
schematically in Fig. 8. 



B. Spin Anisotropy in YBa2Cu3C>6.9 

Our results suggest that a spin anisotropy develops in 
the lower energy (10<-E<30 meV) excitations on entering 
the superconducting state. Nuclear magnetic resonance 
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(NMR) probes the spin fluctuations in the very low fre- 
quency limit and, indeed, the anisotropy of spin-lattice 
relaxation rate (Ti) in YBa 2 Cu30y has been reported to 
show a strong temperature dependence in the supercon- 
ducting state [57, 58]. Various theories have attributed 
this to the combined effect of the NMR form factor and 
a changing x"(q, uj) (See e.g. Ref. 59 and 60). However, 
the present measurements show that there is also an sig- 
nificant intrinsic anisotropy in x^(q, oj) with respect to 
the spin direction which must be considered. It is in- 
teresting to note that Uldry et al. [61] have extracted 
the intrinsic anisotropy from NMR data and concluded 
that the out-of-plane correlations do not change appre- 
ciably on entering the superconducting state, in contrast 
to our results. This may be because NMR measurements 
probe the excitations at much lower frequencies than our 
measurements. 

Anisotropy in the susceptibility ultimately comes from 
the spin-orbit interaction. An exotic case is the super- 
fluid 3 He A-phase [62] , where the susceptibility depends 
on the orientation of the angle of the field to the char- 
acteristic spin vector d. In the case of superconduc- 
tors, dramatic changes in a pre-existing spin anisotropy 
have recently been observed in BaFe1.9Nio.1As2 [36] and 
a small anisotropy at the resonance energy is observed 
in FeSeo.sTeo.s [63]. A possible origin of the emer- 
gence of spin anisotropy in YBa2Cu30g.g may be the 
Dzyaloshinskii-Moriya (DM) interactions between the 
copper spins [64]. The buckled structure of the CuC>2 
planes in ortho-I YBa2Cu306.g (see Fig. 1) means that 
DM interactions of the form D • Si x Sj are allowed 
between neighbouring Cu spins. The presence of such 
terms leads to additional spin anisotropy. This leads to 
a polarization dependence to the spin wave dispersion 
and energy in the antiferromagnetic parent compounds 
La2CuC>4 [65] and YBa2Cu3C>6.2 [66]. In the case of 
YBa2Cu3C>6.2 the anisotropy gaps are ~10 meV [66] and 
the ordered moments lies along the [100] direction [67]. 

The low energy excitations (E <30 meV) we observe 
have their predominant fluctuations within the CuC>2 
planes making the a/b response largest. At higher ener- 
gies, E «40 meV, the excitations are more isotopic. This 
corresponds to all three components of the spin-triplet 
{| ft). I t;>- It), I II)} being excited. 



V. CONCLUSION 

In this work we used inelastic neutron scattering with 
longitudinal polarization analysis to measure the mag- 
netic excitations in the normal and superconducting 
states of near optimally doped YBa2Cu306.g. We have 
unambiguously identified a strong magnetic response in 
the normal state which appears to exist over the 10— 
60 meV range of the present experiment. On entering 
the superconducting state, the out-of-plane magnetic re- 
sponse (Xc)> i s strongly suppressed at lower energies, 
while the response at the magnetic resonance energy and 



above increases. We also find evidence for a new diffuse 
component to the magnetic response in the x'a/b com P°~ 
nent of the susceptibility at low energies 10<£'<30 meV 
which is present in the superconducting state. 
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Appendix A: Sum Rules and the Magnetic Response 

1. Local Susceptibility 

The local susceptibility is a useful way to characterise 
the overall response. It is defined as, 



Sx"(Q,t») d 3 Q 

Jd 3 Q 



(Al) 



where, in general, the integrals are over a volume of re- 
ciprocal space which samples the full Q dependence of 
x"(Q, w). In the case of YBa2Cu3C , 6+ : r this is one unit 
cell in the ah plane and infinity along c. The local sus- 
ceptibility can be split into the two terms of Eq. 5. Thus 
integrating Eq. 5 we have 



x» = x» + x'», 



(A2) 



where 



X » = - / dh / dkx"(h,k,u). (A3) 
1 Jo Jo 

The definition for Xo( w ) us cd here differs by a factor 2 
from earlier work, but allows a direct comparison with 
single layer compounds [11]. 



2. Total Moment Sum Rule 

For an ion with spin only moment, the total squared 
moment is 



= g'f J ,j 3 S(S + l) 

= 3/4 for S = - and g = 2 



(A4) 



The total fluctuating moment observed by INS over a 
given range of energy and momentum can be determined 
from the fluctuation-dissipation theorem and is 



(m 2 ) = (m 2 x + m\ + m 2 ) 



1 — exp(— hw/kT) 



dw. (A5) 
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FIG. 10. (a) Energy-dependent scans of cr£J + crjj at 
Q=(1.5,0.5,1.73) at T = 10,94 K. (b) Difference of scans in 
(a). 



FIG. 9. Wavevector-dependent scans of cyl^ -)- at various 
energies. This combination allows comparison with unpolar- 
ized studies. 



Appendix B: Comparison with Unpolarized Studies 

There are many unpolarized studies of the magnetic 
excitations in YBa 2 Cu306+a; [4, 6, 14, 40, 41]. In this 
section we show that our results are broadly consistent 
with previous results. The main issues that arise in unpo- 
larized studies are: (i) the separation of magnetic signal 
from background and (ii) the separation of magnetic and 
phonon scattering. In the present spin-polarized study 
we may compare to different spin-flip cross-sections to re- 
move the background and the phonon contribution. This 
is demonstrated in Eqs. 3-4. 

The unpolarized inelastic cross section is generally of 
the form 



( d 2 * \ 
\dQdE J 1 - cxp(-hw/kT) 



x"(q,w,T) 



+ N(q,u,T), (Bl) 



where the first term represents the inelastic magnetic re- 
sponse and the second that due to the phonons. A sharp 
magnetic response such as the resonance can be isolated 
through q and u> scans and verified as being magnetic 
through the form factor present in Eq. 1. However, a 
broad or diffuse response is more difficult to distinguish 
from phonons. The phonon response N(q,u,T) usually 



decreases with temperature (Hu> < kT) or remains con- 
stant (Huj 3> kT) due to the Bose factor. Thus a signal 
that increases with decreasing temperature (such as the 
resonance) is likely to be magnetic. If a magnetic signal 
decreases with decreasing temperature] such as the re- 
sponse below about 30 meV in Fig. 3(c)] it is difficult to 
distinguish from phonons using unpolarized neutrons. 

In Figs. 9 and 10, we have reconstructed 'unpolarized' 
scans by adding together the spin-flip and non-spin-flip 
intensities for H || x, ct|J + cr£J:. Our experiment was 
not optimized for this reconstruction because the spin- 
flip channels were counted longer than non-spin-flip, nev- 
ertheless we can make some useful observations. As ex- 
pected, Fig. 9(a) clearly shows the resonance at T = 10 K 
and E = 40 meV in the superconducting state. Note 
there is increased background or phonon scattering at 
larger k in this scan. In the normal state, at T = 94 K, 
it is not possible to identify any magnetic scattering. For 
E = 34 meV [Fig. 9(b)], the scans at both temperatures 
are similar. The data are consistent with a broad mag- 
netic response which changes little between the two tem- 
peratures [see Fig. 5(c)-(d)]. Finally, for E = 26 meV we 
observe a decrease in intensity across much of the scan 
on lowering the temperature. This is consistent with a 
reduction of the magnetic response at this energy [see 
Fig. 5(e)-(f)]. However, the phonon scattering at this en- 
ergy and wavevector is strong [see Fig. 3(a)] thus part 
(about 50%) the reduction observed using unpolarized 
spectroscopy is due to the change of the Bose factor for 
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the phonons. 

Fig. 10 shows energy-dependent scans at the 
Q=(1.5,0.5,1.73) position and a temperature difference 
often used to isolate the resonance (see e.g. [6, 17]). From 



Figs. 3 and 4, we can deduce that about 50% of the ob- 
served change observed with temperature at 26 meV is 
due to phonons. 
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